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Abstract 

We investigate the magnetic properties of nano-holes (NHs) patterned in graphene using first 
principles calculations. We show that superlattices consisting of a periodic array of NHs form a new 
family of 2D crystalline "bulk" magnets whose collective magnetic behavior is governed by inter- 
NH spin-spin interaction. They exhibit long-range magnetic order well above room temperature. 
Furthermore, magnetic semiconductors can be made by doping magnetic NHs into semiconducting 
NH superlattices. Our findings offer a new material system for fundamental studies of spin-spin 
interaction and magnetic ordering in low dimensions, and open up the exciting opportunities of 
making engineered magnetic materials for storage media and spintronics applications. 

PACS numbers: 75.75.+a, 73.21.Cd, 81.05.Uw, 72.80.Rj 
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Magnetic materials have a wide range of applications, such as being used for storage 
media. Magnetism is commonly associated with elements containing localized d or / elec- 
trons, i.e. the itinerant ferromagnetismi^. In contrast, the elements containing diffuse sp 
electrons are intrinsically non-magnetic, but magnetism can be induced in sp-element ma- 
terials extrinsically by defects and impurities. There have been continuing efforts in search- 
ing for new magnetic materials, and much recent interest has been devoted to magnetism 
of carbon-basecr^£&£, especially graphene-based structures^^^^^^^^^^ such as 
graphene nanoribbons^^ 1 ^ 1 ^ and nanoflakes^ 1 ^. Here, we predict a new class of graphene- 
based magnetic nanostructures, the superlattices of graphene nanoholes (GNHs), using first- 
principles calculations. 

Graphene nanoribbonsS&ii^ and nanoflakes^ 1 ^ with zigzag edges have been shown to 
exhibit magnetism. Their magnetization is originated from the localized edge states that give 
rise to a high density of states at the Fermi level rendering a spin-polarization instability 1 . 
Then, for the same reason, if nanoholes (NHs) are made inside a graphene sheet with zigzag 
edges they may also exhibit magnetism. Furthermore, by making an array of NHs, we may 
expect collective "bulk" magnetism because mter-NH spin-spin interactions are introduced 
in addition to the intra-NH spin coupling. This allows us to go beyond the current scope 
limited to the spins within a single nanoribbon or nanoflake. In effect, superlattices consist- 
ing of a periodic array of NH spins form a family of nanostructured magnetic 2D crystals 
with the NH acting like a "super" magnetic atom. 

We have investigated magnetic properties of GNHs, using first-principles pseudopotential 
plane-wave calculations within the spin-polarized generalized gradient approximation^. We 
used a rhombus supercell in the graphene plane with the cell size ranging from 14xl4A to 
41x41A and a vacuum layer of ~10A. We used a 2x2x1 k-point mesh for Brillouin zone 
sampling and a plane wave cutoff of 22.1 Rd. The systems contain up to a maximum of 530 
atoms. The dangling bonds on the edge atoms are saturated with hydrogen. The system is 
relaxed until the force on each atom is minimized to less than 0.01 eV/A. 

Considering first a single zigzag NH by examining the intra-NH spin-spin interaction, we 
found that individual NH can be viewed as an "inverse structure" of nanoflake or nanoribbon, 
like an anti-flake or anti-ribbon, with similar spin behavior. We determine the ground-state 
magnetism of three typical NH shapes: triangular (FigfjJ,), rhombus (FigJTb) and hexagonal 
(Fig {lb), by comparing the relative stability of ferromagnetic (FM), antiferromagnetic (AF) 
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and paramagnetic (PM) configuration as a function of NH size. Our calculations show 
that the ground state is FM for triangular NHs, but AF for rhombus and hexagonal NHs, 
and their corresponding spin densities are shown in Figflh, db andTb, respectively. The 
magnetic moments are highly concentrated on the edges and decay quickly away from the 
edge, as shown in FigflU. Similar decaying behavior has been seen in nanoribbons&ii and 
nanoflakes^ 1 ^. The edge moment increases with increasing NH size (see inset of Fig{TJi). 

The triangular NHs have a metastable ferrimagnetic state with two edges having one 
spin and the other edge having the opposite spin (see supplementary Fig{T]). For a 4-atom 
triangular NH (FigHh), the FM state is 52 meV lower in energy than the ferrimagnetic state, 
and the latter is 13 meV lower than the PM state. For a 32-atom rhombus NH (FigfTb), 
the AF state is 89.2 meV lower than the PM state; for a 54-atom hexagonal NH (Fig{Tb), it 
is 164.4 meV. The energy difference increases with increasing NH size. The triangular NHs 
favor FM at all sizes, whereas rhombus and hexagonal NHs only become AF when the edge 
contains more than five atoms, i.e. they are PM if the NH is too small. So, the triangular 
NHs have a stronger tendency toward magnetization. 

The magnetic ordering within a single NH is consistent with both the theorem of itinerant 
magnetism in a bipartite lattice^ and the topological frustration model of the 7r-bonds 17 
counting the unpaired spins in the nonbonding states^ 1 ^. For a system like graphene con- 
sisting of two atomic sublattices, each sublattice assumes one spin and the total spin S of the 
ground state equals \\N B — Na\ where N B (N^) is the number of atoms on B (A) sublattice. 
Because of the honeycomb lattice symmetry, atoms on the same zigzag edge belong to the 
same sublattice; while atoms on two different zigzag edges belong to the same sublattice if 
the two edges are at an angle of 0° or 60°, but different sublattices if at an angle of 120° or 
180°. Consequently, the triangular NH are FM, because all three edges are 60° to each other 
on the same sublattice; the rhombus and hexagonal NHs are AF, because one-half the edges 
are on the A-sublattice and another half on the B-sublattice as the two types of edges are 
120° to each other. This same argument can be applied to nanoribbons 9 and nanoflakes^ 1 ^. 

Next, we consider GNH superlattices (a periodic array of NHs) by examining the inter- 
NH spin-spin interaction. In principle, one can generate four out of five possible 2D Bravais 
lattices of NHs (see supplementary Figj2]). Here, we focus on the honeycomb superlattices of 
triangular NHs (Figj2^i and [2b), in which each NH possesses a net moment acting effectively 
as "one" spin. The superlattice contains two sublattices of NHs, superimposed on the 
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FIG. 1: The ground-state magnetic configurations of different shapes of NHs. (a) FM triangular 
NH; (b) AF rhombus NH; (c) AF hexagonal NH. In (a-c), white and red balls indicate the up- and 
down-spin density isosurface at 0.02e/A3, respectively; blue and white sticks represent C-C and 
C-H bonds respectively, (d) The average local magnetic moment (hb) per atom in the triangular 
NH (FigJTh.) as a function of distance moving away from the center of NH, measured in atomic 
shells with the edge atoms as the first shell. The inset shows /ib on the edge vs. NH size (I). 

background of graphene containing two sublattices of atoms. We realize that the NHs on 
the same sublattice will be FM-coupled because their corresponding edges are at 0° to each 
other so that their edge atoms are on the same atomic sublattice. On the other hand, the NHs 
on different sublattices will be FM-coupled if they are in a parallel configuration (FigfS^) but 
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AF-coupled if they are in an antiparallel configuration (Fig|2b) when their corresponding 
edges are at 180° to each other so that their edge atoms are on different atomic sublattices. 
These have indeed been confirmed by our first-principles calculations. Independent of NH 
size and supercell dimension, the FM state is favored for parallel configurations but the AF 
state is favored for antiparallel configurations. In both cases, the spin-polarization splits the 
edge states opening a gap at the Fermi energy- 1 ^. The total spin S in one unit cell equals 
to l\N B -N A \; it increases linearly in the FM parallel configuration but remains zero in the 
AF antiparallel configuration with increasing NH size. 

The collective magnetic behavior of a GNH superlattice depends on inter-NH spin-spin 
interaction. Particularly, there exists super exchange interaction between the NH spins, in 
addition to the spin coupling defined by the underlying bipartite lattice (i.e., the relative 
angle between the zigzag edges of GNHs). In FigfSb, we plot AE pc = E(FM) — E(PM) for 
the FM parallel configuration and AE ac = E(AF) — E(PM) for the AF antiparallel con- 
figuration as a function of cell dimension (L), i.e., the NH-NH separation. |A.E pc | increases 
while |Ai? ac | decreases with decreasing L. This indicates that as the NHs getting closer, the 
FM state becomes relatively more stable than AF state in both configurations, i.e. the FM 
coupling is favored by the super exchange interaction. Also plotted in FigfSb are magnetic 
moments on the NH edges, which are found to increase in the FM but decrease in the AF 
configuration with decreasing L. This again reflects that the edge magnetization on the 
neighboring NHs is enhanced with the same spin when they are FM coupled but suppressed 
with the opposite spin when they are AF coupled by the super exchange interaction. 

The above results show that long-range ferromagnetic ordering can be created by employ- 
ing the parallel configuration of triangular NHs in different lattice symmetries, as illustrated 
in supplementary Fig. [21 The next important question is what Curie temperature (T c ) they 
can have. We have estimated T c using the mean-field theory of Heisenberg model^i^ 1 ^, 



Where A is the energy cost to flip one "NH spin" in the FM lattice, which have been 
calculated directly from first principles for the honeycomb lattices (Fig 12k,). Figure [2jd shows 
that T c increases from 169 K to 1388 K when NH size (I) increases from 0.738 to 1.476 
nm with cell dimension (L) fixed at 2.982 nm, and decreases from 586 K to 169 K when L 
increases from 1.704 nm to 2.982 nm with / fixed at 0.738 nm. These trends are expected since 
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FIG. 2: (a) Ground-state spin configurations in a FM honeycomb NH superlattice. (b) Same as 
(a) in an AF superlattice. All the symbols and notations for bonds and spin densities are the same 
as FigHJ Yellow dashed lines mark the primitive cell, (c) AE pc = E(FM) - E(PM) of the FM 
superlattice (red solid squares) and AE ac = E(AF) — E(PM) of the AF superlattice (blue solid 
triangles) versus cell dimension (L); Edge magnetic moments, in the FM lattice with fixed hole 
size (1=1.476 nm) (red open squares) and in the AF lattice (Z=0.738 nm) (blue open triangles) 
versus L. (d) Curie temperature of the FM superlattice as a function of NH size (I) and L. 

magnetization is stronger for larger NH size and higher NH density. Limited by computation 
time, some of our cell dimensions are possibly unrealisticaliy too small (NH density too high), 
which gives rise to a very high T c . Still, it is important to point out that our calculations 
suggest that it is possible to make FM GNH superlattices with T c above room temperature 
by using a NH size of 50 nm and a density of 10 -4 nm -2 , achievable by today's lithographic 
patterning technology. We note that a recent experiment^ 4 - has shown a T c > 350 K in FM 
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FIG. 3: Illustration of DMS made from GNH superlattice. (a) A semiconductor GNH hexagonal 
lattice (L = 8y/3a,a = 2.46A is the lattice constant of graphene.) consisting of an array of rhombus 
armchair NHs (L = 8y3a). (b) GGA Band structure of (a), the inset shows the Brillouin zone, 
(c) TB band gap of (a)-type structures as a function of NH size (I) and cell dimension (L). (d) 
Magnetic semiconductor made by doped (a) with triangular zigzag NHs. All the symbols and 
notations for bonds and spin densities are the same as Fig(TJ 

graphite made by proton bombardment. 

Since graphene-based nanostructures hold great promise for future electronics^ 1 ^ 1 ^ 1 ^, 
our discovery of GNH magnetism offers the exciting prospect of combining magnetic and 
semiconducting behavior in one material system. Here, we demonstrate the possibility of 
making diluted magnetic semiconductors (DMS) by exploiting GNHs with two different 
kinds of edges. Similar to superlattices of zigzag NHs, we can create superlattices of armchair 
NHs, which constitute a class of 2D semiconductors. Figure [3]d shows the band structure 
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of a superlattice of rhombus armchair NHs (FigO^i) having a direct band gap of 0.43 eV, 
as obtained from first-principles calculations. Figure &p shows the band gap as a function 
of NH size (/) and cell dimension (L), from tight-binding calculations 3 -. The gap increases 
with increasing I but decreases with increasing L. 

DMS can be made by adding triangular zigzag NHs into the semiconductor superlat- 
tice, as illustrated in Fig. [3H. To ensure the ferromagnetism, all triangular NHs must be 
parallel with each other acting like magnetic dopants. Usually DMS are synthesized by mix- 
ing two different materials, typically III-V semiconductors and transition-metal magnetic 
element s^ 1 ^. The main challenge is to increase the magnetic dopant concentration in order 
to raise the Curie temperature, because the two types of materials are usually not miscible. 
Here, we introduce an "all-carbon" DMS, in which combined semiconductor and magnetic 
behavior are achieved by structural manipulation. Consequently, room-temperature DMS 
are possible because the dopant concentration can be increased without the miscibility prob- 
lem. One might also consider doping other magnetic elements into the semiconducting GNH 
superlattice. 

It is very exciting to consider making engineered magnetic materials with NHs for various 
applications. For example, it is possible to directly pattern NHs into engineered magnetic 
storage media (see Fig. H]). The ground state of rhombus NHs is AF (Fig. [lb and FigJH 
lower-left inset) and their first excited state is FM (Fig. HI up-right inset) when the NH size 
is larger than 14. 6 A according to our calculation. Taking each NH as one bit, we can assign 
the ground state with "S=0" and the excited state with "S=N" to represent the '0' and 
'1', respectively. The switching between '0' to T' can be done by applying a local magnetic 
field or energy pulse to convert between the ground and the excited state. Using a NH size 
of 50 nm and a density of 10 -4 nm" 2 , a storage density about 0.1 terabit per square inch 
would be achieved, much higher than the current density in use. One interesting topic of 
future study is the magnetocrystalline anisotropy around individual NHs, which must be 
larger than ksT for the proposed storage media to work. 

We thank DOE-NERSC and Center for High Performance Computing (CHPC) at the 
University of Utah for providing the computing resources. This work was supported by 
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FIG. 4: Schematic Illustration of a magnetic storage medium consisting of a patterned array of 
rhombus GNHs. The insets show the detailed structure of "0" and "1" bit, represented by the 
ground-state AF configuration (S=0) and the excited FM configuration (S=N), respectively. Red 
and green balls show the spin-up and spin-down density at an isosurface value of 0.02e/A 3 . 
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Supplementary Fig.l The spin-density plot of the ferrimagnetic configuration of a 4-atom 
triangular NH. White and red balls indicate the up- and down-spin density isosurface at 0.02e/A 3 
respectively; blue and white sticks represent C-C and C-H bonds respectively. 
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Supplementary Fig. 2 Schematic illustration of four possible types of Bravais lattice of GNHs that 
can be patterned in graphene. Solid arrows and lines mark the primitive cells, (a) hexagonal 
lattice; (b) rectangular lattice; (c) centered rectangular lattice; the dashed lines mark the 
conventional cell; (d) oblique lattice. Note that the square lattice is not possible. 
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